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Abstract 

The zero-temperature phase diagram of binary mixtures of particles interacting via a screened 
Coulomb pair potential is calculated as a function of composition and charge ratio. The potential 
energy obtained by a Lekner summation is minimized among a variety of candidate two-dimensional 
crystals. A wealth of different stable crystal structures is identified including A, B, AB2, A2B ^ AB4 
structures [A (B) particles correspond to large (small) charge.] Their elementary cells consist of 
triangular, square or rhombic lattices of the A particles with a basis comprising various structures 
of A and B particles. For small charge asymmetry there are no intermediate crystals besides 
the pure A and B triangular crystals. The predicted structures are detectable in experiments on 
confined mixtures of charged colloids or dusty plasma sheets. 
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I. INTRODUCTION 



Two-component mixtures in general exhibit much richer crystalhzation phenomena and 
polymorphism than their one-component counterparts 1{| as witnessed by a hiige variety of 
possible stable binary crystals, e.g. for binary hard sphere systems [2|, y, y, l5j. How the 
whole crystal phase behavior in mixtures depends on the interparticle interactions is far 



from being understood even in equilibrium 



aepe 

m 



This is true also in two spatial dimensions 



where the number of Bravais lattices is smaller than in three dimensions. Binary mixtures in 
two dimensions have been studied for hard disks [s] and a complex diagram of close packing 
was obtained as a function of their diameter ratio. More recently, a two-dimensional binary 
mixture with soft interactions was considered [9| , namely that for parallel dipoles where the 
pair potential scales with the inverse cube of the interparticle separation. A variant of this 

granular matter 



model has been considered in Ref. [10|. Such systems can be realized in 



ll| and in magnetic colloidal suspensions confined to a air-water interface 



12] ■ Again, as a 



function of the ratio of dipole moments of the two species, a complex stability phase diagram 
of stable binary crystals was obtained that qualitatively differs from the hard disk case j^. 
In particular for low asymmetries, the hard disk system shows a complete separation into 
pure A and B triangular crystals [Sj while the soft dipolar systems possesses two stable 
mixed crystals as well with stoechiometric ratio A2B and AB2 9|. These differences show 
that the topology of the phase diagrams depend on details of the interactions and there is 
certainly a need to understand this dependence in more detail. 

In this paper, we consider a two-dimensional binary system of Yukawa particles, i.e. 
the pair interaction potential V{r) between the particles is a screened Coulomb interaction 
oc exp(fi;r)/r where k is the screening constant (or the inverse screening length). This 
potential interpolates between the case of hard disks (as obtained in the limit of high k) 
and the unscreened Coulomb case (as obtained for k, = 0). The latter limit, V{r) oc 1/r is 
even softer than the dipolar case where V{r) oc The two components are defined by 

two different charges, i.e. different prefactors in front of the Yukawa interaction. In previous 
works, such a classical binary mixture with Yukawa interactions in three-dimensions has been 
used as a model to study mixing rules 



15l . Il6l . [iTJ, dynamical correlations |l8l . Il9l | and transport properties 



3|, effective forces [14(], fluid-fluid phase separation 

2(jl. Likewise the 



pure (one-component) Yukawa system was also studied in two-spatial dimensions for fluid 
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structure 
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2J], dynamics 
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271 . |28| and transport properties 29|. Binary 



mixtures of Yukawa particles in two dimensions have also been studied for fluid structure 



30|, adsorption [3l|, interfaces j32|] and transport [33|. However, the crystallization issue 



was only addressed in one-component Yukawa systems (for a recent work, see e.g. 3J]) but 
never in binary mixtures. 

The Yukawa potential is realized in charged colloidal suspensions jsS^ as well as in dusty 
plasmas [36], both for one component systems and mixtures. In fact, highly charged col- 
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39| 



loidal suspensions can be confined between highly charged parallel glass plates [37 
which restricts their motion practically to two dimensions. The interactions between these 
macroions are screened due to the presence of the microscopic microions and additional salt 



ions. As in three dimensions, the Debye-Hiickel screened Coulomb interaction is a reason- 
able model for confined charged colloids 40, 4l|. Crystallization of binary charged colloids 
has been studied experimentally in the bulk. However, a monolayer of a confined binary 
mixture of charged colloids has not yet been realized although this is in principle possible as 
has been shown for sterically-stabilized 421] and magnetic colloids 43|. On the other hand, 
sheets of highly charged dust particles in plasmas (so-called complex plasmas) can also be 
confined to two dimensions, e.g. by levitating electric fields. The interaction between the 
dust particles is again screened such that a Yukawa model is appropriate [Sg, 14^, |45|. Highly 



charged microspheres suspended in a plasma settled in a horizontal monolayer were studied 
experimentally and compared to a two-dimensional Yukawa model 0, 47, 48] . There is no 
principle problem in studying binary mixtures of dust particles but a concrete realization in 
an experiments still has to be performed as well. 

Apart from its important realizations, our major motivation for our studies is to un- 
derstand the interplay between the interparticle interaction and the stability of different 
two-dimensional crystal lattices. A control of colloidal composite lattices may lead to new 
photonic crystals 49 ] to molecular-sieves |50!] and to micro- and nano- filters with desired 
)orosity [5l|- The electric properties of a nanocrystal depend on its superlattice structure 
52l |. For these type of applications, it is crucial to understand the various stable lattice 
types in binary mixtures. 

For the two-component two-dimensional Yukawa mixture, we obtain the full phase di- 
agram at zero-temperature as a function of the charge asymmetry using lattice sums. 
As a result, we find a variety of different stable composite lattices. They include 
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A, B, AB2, A2B, AB4 structures. Their elementary cells consist of (equilateral) triangular, 
square and rhombic lattices of the big particles. These are highly decorated by a basis in- 
volving either A particles alone or both B and A particles. The topology of the resulting 
phase diagram differs qualitatively from that of hard disk mixtures ja] and dipoles . 

The paper is organized as follows: In Sec. II the model is described and possible candidate 
structures for crystal lattices in two dimensions are proposed. Results for the phase diagrams 
are presented in Sec. III. We conclude finally in Sec. IV. 



II. MODEL 



The model systems used in our study are binary mixtures of (repulsive) charged particles 
made up of two species denoted as A and B. Each component A and B is characterized 
by its charge valency Za and Zb, respectively. These constitutive particles are confined to 
a two-dimensional plane, and interact via the Yukawa pair potential. Introducing the ratio 
Z = Zb/Za the pair interaction potentials between two A particles, a A- and 5-particles, 
and two -B-particle at distance r are 

VAA{r)= KVoip{r), VAsir) = KVoZLp(r), 

VBB{r) = KVoZ'ifir), (1) 

respectively. The dimensionless function <f{r) is given by 

^(r) = ^^3^. (2) 

K,r 

where the energy amplitude VqK sets the energy scale. 

Our goal is to determine the stable crystalline structures adopted by the system at zero 
temperature. We consider a parallelogram as a primitive cell which contains ua A-particles 
and ub -B-particles. This cell can be described geometrically by the two lattice vectors 
a = a(l, 0) and b = a7(cos 6, sin 6), where 6 is the angle between a and b and 7 is the aspect 
ratio (7 = |b|/|a|). The position of a particle i (of species A) and that of a particle j (of 
species B) in the parallelogram is specified by the vectors rf = {xf, yf) and = {x^, y^), 
respectively. The total internal energy (per primitive cell) U has the form 
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1 



4 - vj + R|) 



J=A,B i,j=l R 



(3) 



i=l j = l R 



where R = ka + Ih with /c and / being integers. The sums over R in Eq. [3] run over all 
lattice cells where the prime indicates that for R = the terms with i = j are to be omitted. 
In order to handle efficiently the long-range nature of the Yukawa interaction at moderate 
screening strength, we employed a Lekner-summation (see Appendix A). 

We choose to work at prescribed pressure p and zero temperature (T = 0). Hence, 
the corresponding thermodynamic potential is the Gibbs free energy G. Additionally, we 
consider interacting particles at composition X := ^^/(n^ + ns), so that the (intensive) 
Gibbs free energy g per particle reads: g = g{p,Z,X) = G/{nA + tib)- At vanishing 
temperature, g is related to the internal energy per particle u = U/{nA + ^b) through 
g = u + p/ p, where the pressure p is given hj p = p'^{du/dp), and p = {ha + nB)/\sL x b| 
is the total particle density. The Gibbs free energy per particle g has been minimized with 
respect to 7, 6 and the position of particles of species A and B within the primitive cell. 
In order to decrease the complexity of the energy landscape, we have limited the number of 
variables and considered the following candidates for our binary mixtures: A4B, A3B, A2B, 
A^Bi, A3B2, AB, A2B2, A^B^i, A2B3, AB2, ^2^4, AB3, ABi and ABf^. For the ABq and 
A-^B^ case we have only considered a triangular lattice formed by the A particles. 

III. RESULTS 

A. Phase diagram 

The ultimate phase diagrams in the (Z, X) plane has been obtained by employing the 
Maxwell construction. We recall here that the both dimensionless quantities, namely the 
charge ratio Z as well as the composition X, can vary between zero and unity. A low 
charge ratio (i.e., Z is close to zero) indicates a strong charge asymmetry, whereas a high 
charge ration (i.e., Z is close to unity) represents a large charge symmetry or equivalently 
a weak charge asymmetry. Given the fact that the phase behavior is getting increasingly 
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T(A)A^B^ Rh(A)AS2 Rh(A)AS 




FIG. 1: The stable binary crystal structures and their primitive cells. The red (green) discs 
correspond to A (B) particles. 
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TABLE I: The stable phases with their Bravais lattice and their basis. 



Phase Bravais lattice [basis] 

T{A) Triangular for A [one A particle] 

T{B) Triangular for B [one B particle] 

S{AB) Square for A and B together [one A and one B particles] 

S{A)Bn Square for A [one A and n B particles] 

'Rh{A)AraBn Rhombic for A [(m + 1) A and n B particles] 

T{A)AmBn Triangular for A [(m -\- 1) A and n B particles] 



complicated upon lowering Z, involving a huge basket of candidates, we only present results 
starting from Z = 0.2. Furthermore, in contrast to situations where the pair potential can 
be described as a power low of the separation distance (as it was the case in our previous 
work on dipolar mixtures [9]), the phase diagram becomes pressure dependent for Yukawa 
systems. To capture this feature, we present results at three well distinct pressures, namely 
p* = p/[VoK^) = 0.01, 1 and 100. An overview of the resulting stable crystalline phases can 
be found in Fig. [H The corresponding nomenclature of the phase labeling is explained in 
Table [B The phase diagrams in the {Z, X) plane for the three reduced pressures p* = 0.01, 1 
and 100 are depicted in Fig. [2]^a), Fig. Mjo), and Fig. MyC), respectively. Note that upon 
increasing p* at prescribed Z and X one decreases the density. 

Let us first focus our discussion on the apparently simple phase behavior reported at 
weak charge asymmetry (here roughly Z > 0.5, see Fig. [2]). Thereby, the system phase 
separates into a pure A-triangular crystalline phase and a B-one (see also Fig. [1]). This 
triangular structure obviously corresponds to the single-component ground-state. Having 
in mind that the same phase behavior is reported for hard disks binary mixtures at small 
size asymmetry [sj, it is meaningful to equally expect a phase separation for moderate or 
sufficiently large reduced screening strength k* = ^^1 \fp- For Z = 1, we have k* ^ 3.0, 1.2, 0.4 
for p* = 0.01, 1, 100, respectively, so that the phase separation is certainly to be expected 
for moderate pressures (here p* = 0.01 and possibly p* = 1) when referring to the hard disk 
hmit fl. 

What is now less obvious, still in the regime of weak charge asymmetry, is the phase 
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FIG. 2: The phase diagram in the (Z, X) plane of charge asymmetry and composition at T = for 
a effective pressure (a) p* = 0.01, (b) p* = 1, (c) p* = 100. The symbol # (*) denote continuous 
(discontinuous) transitions. 

separation reported in Fig. Ws) for p* = 100, where k*^ ~ 0.4 for large charge symmetry. 
Recently, we have shown for dipolar binary mixtures |9;], whose pair potential is governed 
by that, at weak dipolar asymmetry (the analogous quantity to the charge ratio in 

our present study), the stable mixtures A2B and AB2 (who are globally triangular) set in. 
This phase behavior contrasts therefore strongly with that presently reported for Yukawa 
mixtures, see Fig. [2](c). Given the fact that at weak screening the Yukawa pair potential 
is well approximated by a 1/r dependence, which is even softer than 1/r^, it is legitimate 
to expect stable mixtures in the regime of weak screening and charge asymmetry. In order 
to check this idea we have performed additional calculations at p* = 10^° with Z = 0.99 
leading to reduced screening strengths of the order of 10~^. Those values for k* turn out 
to be still too large to recover the phase behavior found at 1/r^-pair interactions j^. The 
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consideration of even much smaller screening strengths (say roughly of the order of 10~^) are 
numerically not tractable within reasonable CPU time. Unfortunately, the implementation 
of a direct Lekner and/or Ewald sum for the 1/r-pair interactions is delicate at prescribed 
pressure, since the lack of electroneutrality involves the presence of an artificial homogeneous 
neutralizing background which is thermodynamically only consistent at prescribed density 



531]. Consequently, although we have a strong intuition about the stability of mixtures at 



weak charge asymmetry and screening, we can not prove it here on computational basis. 

We now briefly address the more complicated phase behavior reported at strong charge 
asymmetry, see Fig. [2] with Z < 0.5. As a clear general trend, it is found that the number 
of stable phases increases with growing pressure. This feature is in agreement with the idea 
that mixing is favored upon softening the pair potential. 

A common and remarkable feature in this regime of strong charge asymmetry (see Fig. [2]) 
is the imposing stability of the composition X = 1/2. This feature was also reported for dipo- 
lar mixtures [9|. More specifically, the following cascade S{AB) T(A)A2B-i Iih(A)AB2 
is found upon increasing Z, see Fig. |2]and Fig. [T]for the corresponding structures. Thereby, 
the transition S{AB) -> T{A)A2B3 is discontinuous whereas T{A)A2B3 -> Rh{A)AB2 is 
continuous, see Fig. [2l Note that, for p* = 0.01 shown in Fig. [2]^a), the stability of the 
square phase S{AB) occurs for values of Z smaller than 0.2 that are not shown here. 



B. Thermodynamical properties 

1. Constant pressure 

In this part, we investigate some thermodynamic properties such as the reduced density p* 
or the reduced Gibbs free energy g* = g/iVoK,), as obtained prior the Maxwell construction. 
Although the pressure considered here is fixed at p* = 100, very similar results are obtained 
for the two other pressures. 

The reduced density p* as a function of the charge ratio Z at different compositions X 
is sketched in Fig. [31 At given composition X, the density decreases monotonically with Z, 
see Fig. |3l This effect can be simply explained as follows: Upon increasing Z the repulsive 
A — B and B — B pair interactions increase accordingly, so that to keep the pressure fixed the 
system has to decrease its density. Moreover, at prescribed charge ratio. Fig. Vindicates that 
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0.2 0.4 0.6 0.8 1 

z 

FIG. 3: Reduced density p* (prior the Maxwell construction) as a function of the charge ratio Z 
for various compositions X at prescribed reduced pressure p* = 100. The arrow indicates growing 
X. 

the density increases with the composition. This feature can also be explained with simple 
physics: Upon enlarging the composition X, the proportion of weakly charged 5-particles 
increases accordingly, so that to keep the pressure constant the system has to increase its 
density. 

The reduced Gibbs free energy g* as a function of the charge ratio Z at different com- 
positions X is sketched in Fig. HI Recalling that g* = u* +p*/p* [with u* = u/iVoK)], the 
behavior of g* exhibi ' ~ D seating the just 
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FIG. 4: Same as Fig. E]but for g*. 
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FIG. 5: Reduced density p* (prior the Maxwell construction) as a function of the charge ratio Z 
for various reduced pressures p* at prescribed composition X = 1/2. 

explained behavior of p* described in Fig. [3l Since the reduced internal energy u* decreases 
with growing X and therefore with growing p* at prescribed Z, it is clear that g* decreases 
with growing X at given Z, as seen in Fig. |H Besides, at given composition X, Fig. H] 
shows that g* increase with growing Z, as expected. 

2. Constant composition 

We now analyze p* and at X = 1/2 as a function of Z for different values of p* . As 
far as the behavior of p* is concerned, the new information provided by Fig. [5] is that p* 
increases with growing pressure at given charge ratio, as it should be. The same qualitative 
feature is also observed for g* in Fig. [HI A closer inspection of Fig. O and Fig. [H] suggests 
that p* and g* increase rather slowly with p* (at given Z*). 

IV. CONCLUDING REMARKS 

In conclusion we have determined the ground-state (i.e. zero-temperature) phase diagram 
for a two-component Yukawa monolayer at various pressure for arbitrary compositions and 
a broad range of charge asymmetries. Among a big number of candidate phases, a wealth of 
different composite lattices has been found to be stable. The larger the charge asymmetry, 
the more complex is the phase diagram. At low asymmetry the system shows demixing 
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FIG. 6: Same as Fig. [5] but for g* . 



into pure A and B crystals similar to hard disks but different from the soft inverse cube 
interaction valid for dipoles. The results are in principle detectable in binary mixtures of 
charged colloids confined between two charged plates or levitated dusty plasma sheets. 

It would be interesting to study the effect of finite temperature. We expect that the 
topology of the phase diagram does not change upon gently increasing the temperature 



though this could change close to melting. When cooling a two-component fluid down, g. 
formation in the binary systems at finite temperature may be a fascinating topic as well 



ass 
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to be studied in the future. In fact, it has been speculated that the underlying crystallization 



into the stable crystal lattices may control vitrification 
relevant for the structure of glasses as well. 



55| and therefore are findings are 
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APPENDIX A: LEKNER SUMS FOR YUKAWA INTERACTIONS IN TWO DI- 
MENSIONAL SYSTEMS 



We consider a primitive cell in the shape of a parallelogram, which contains a set of 
n = riA + ub particles interacting via Yukawa potentials. The parallelogram repeated in 
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the xy plane gives a 2- dimensional lattice, and can be described by two lattice vectors 
a = (a^,0) and b = {bx,by). In the parallelogram, the position of a charge valency Zi is 
defined by = {xi,yi). 

The total interaction energy per cell is given by 



U _ 1 
Vo~2 



i=i 



i=l 



with 



where 



>r;^ exp(-K|r + R|) exp(-/s;|R| 

R ' ' R^O ' ' 



(Al) 



(A2) 



|r + R| = \J [x + aj + hxTTiY + (?/ + hymY and |R| = \J {aj + hxmY + {byiny. 



Here R = /a+mb with / and m being integers. The slowly convergent sums over lattice sites 
(Eq. IA2p can not be efficiently used in a numerical calculation, so that we will transform 



them into rapidly convergent forms using a Lekner Method 
following integral representation 

.2 



56, 



57|. With the help of the 



exp( 


— K 


r + R|) 




r + R 





- , ^exp(-— - |r + R| t), 
TT Jo y/t 4t 



(A3) 



we obtain 



1 f°° dt { \ ( X 



Now, to get further, we apply a 1-dimensional Poisson summation 



oo oo . 

^ exp [-(a + /5/)'t] = ^ 5^ exp 

l=—oo '"^ k=—oo ^ 



a \ Ti 



(A4) 



(A5) 



which provides 

+00 



a exp 

i=— oo 



-(— + / + m— )2a^t) 

dry diry 



1 /vr 



1 + 2E 



COS 



A:=l 



27r/c ( h m— 



exp (— TT^/c^/a^i) 
(A6) 



14 



Inserting Eq. (1A6[) into Eq. (1A4[) yields: 
^ 1 >^ r cit 



r?i=— oo 

+00 +00 

' fc=l m=— cxD 

-exp 



cos 



27rA; ( h m — 



2U2 



A7f'k'\ 1 



— -(?/ + mby) t 



Now, taking into account the following relation 



dt 
1 



exp 
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' 4t 



CH = 2Ko{BC) 



(A7) 



(A8) 



where Kq is the zeroth order modified Bessel function of the second kind, The final expression 
for $(r) reads: 



+00 



<l>(r) = ^ 5Z Ko{K\y + mby\) 

lO-x I 

4 



m=— 00 

00 +00 

fc=l m,=— 00 



COS 



, X bx 
27r { h m — 



xKo 

and the "self" contribution $0 

00 



I y + mby \\ 



47r2A;2 



for 2/ 7^ 



(A9) 







^ V — V 

I — ;2^Ko{Kmby) 

l*^^'' m=l 

— \ 271""''^ — I -^0 

' k=l m=l ^ ' 



2 



mb^ \ K? + 



In [1 — exp(— fta^:)] 



(AlO) 



In the limit of a rectangular based cell, i.e setting b^ = 0, one obtains the formulas for the 
cross and self-energies that are identical to those derived in 57|] with z = 0. 
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